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Тестування якост! обчислення оцнки 
перетворення Фур’е фтштних функшй 


У статт! проведено тестування якост! алгоритм обчислення перетворення Фур’е на класах функций, що 
задовольняють умову Лшшиця, та в1дпов1дних 1м 1нтерполящйних клас1в та показано шдвищення якост! 
обчислень шляхом занурення задач! в б1лыш вузький клас. Були отриман! апостерорн! ощнки точност! для 
розглянутих тестових задач. Була шдтверджена шляхом тестування яюсть квадратурних формул та 
апр!орних ощнок 1х точност1. У статт! наводяться фрагменти таблиць з результатами тест. 


Вйдомо, що розробка великих комплексов прикладних програм для розв’язування 
складних математичних задач — досить складна проблема. Шсля забезпечення загально! 
працездатност! системи виникае задача дослдження цих систем з метою виявлення 1х 
реальних показникв якост! та можливостей застосування в процес! оптимзаци обчис- 
лень. Важливсть рол! тестування алгоритмв-програм загалом та чисельного експеримен- 
ту зокрема в цих процесах була обтрунтована в ряд! праць [1], [2], крм того, комп’ютерна 
технолопя розв’язання задач з шдвищеною якстю [3] базуеться саме на результатах 
тестування алгоритмв-програм. 

Перетворення Фур’е широко застосовуеться в цифровйй обробщ сигналйв та зобра- 
жень. Тому задача виявлення апостерорних ощнок, як одного з резерыв оптимпзацй, для 
алгоритмв обчислення перетворення Фур’е на р1зних класах задач шляхом тестування 
досить актуальна. 

Вдомо, що одним з резерыв оптимзаци е занурення задач! в б1лльш вузький клас з 
максимальним використанням апр!орно! 1нформацй. Мета дано! стати — практична 
тлюстрацщ використання цього резерву оптим1зацй обчислень шляхом тестування. 


Постановка задач! тестування 


Дамо формальне означення оптимального за точнистю алгоритму [4]. 
Нехай замисть елемента /, що належить до множини вхдних даних задач! Р, маемо 


ск!нченного вимру вектор Гу, = (я д —| ‚ а = характеризуе точнсть фактичних значень 
вхдних даних [у =(1,...йу): В 1 |< &, /=1, М. Позначимо К = АГ, , -— наближений 
розв’язок задач1 Р (2 ‚ К) ‚ де КЕ - точний розв’язок задачт, А — алгоритм, застосо- 


ваний до [,,. Тод! абсолютною похибкою наближеного розв’язку задач! Р буде 


(А, №,=,Р)= ря (®. В ) Введемо характеристики 
Г (А.М, =,Р)) = зирзир! (А, №,=,Р); 
РеЕР Ке\ 
(М, М, =,Р) = шШЕУ (4, М,=,Р); 
Г=И(М,=,Р) =шЕТ (4, М,=,Р); 
А 


(м№,Р) =У (М№,0,Р), 
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де Р- вдповдно заданий клас задач Р, М — множина алгоритмив А розв’язування задач, а 
(М ‚=,Р) — нижня границя по вом алгоритмам А. Алгоритм А, на якому досягаеться 


7 (№ ‚5, Р) ‚ назвемо оптимальним за точнистю алгоритмом. 
Якщо для алгоритму А 
ИА, М, =,Р}<ИМ,=,}+1, п>0, 


то А назвемо оптимальним за точнстю з точнстю до 7. Якщо и = о[И| або ОГИ|, то А 
назвемо вдпов1дно асимптотично оптимальним або оптимальним за порядком. 
Потрибно обчислити 1нтеграл в1д швидкоосцилюючо! функцй виду 
Ь 


= | У(х)а(ех)ах, (1) 


а 


2л 
де Г (х) РЁ (Е- деякий клас функц!И), @ — довйльне дйсне число, @| 2 6—2 1, 1нформа- 


щя про Хх) задаеться в М вузлових точках 3 [а, Ь]. 
В робот! вибран! класи С, (Лшшиця) та вдповдний йому 1нтерполящйний клас 
С, х › що езвуженням попереднього: 
С, — клас функщй Хх), визначених на [а, 6] 1 таких, що задовольняють умову: 
и) 6 < [4,5]. (2) 


С, х — клас функций Г(х)еС, тзаданий фиксованими значеннями у вузлах фиксовано! 


сики х, =0,М-1. 


Оптимальн! за точнистю квадратурн! формули та апр1орн! 
ощнки похибки 


ВИдповдно до теорем, доведених в [5], для класу Липпица оптимальною за точнстю 
буде квадратурна формула 


— +112 


Ул [за(охах, де хану = + АА (3) 
1=0 . 


Х-1/2 


1 ощнка похибки 


Г У 
Исьо т Уз [7 [бул | +76) (4) 
[2] =0 2 
де 
А @Ах м Ах м1 
Р(®)= - С 2 }=| 5 ] (1 хм-1)с0$(6). (5) 
Для класу Сьм оптимальна квадратурна формула матиме вигляд: 
М1 
И | Г“ зна (ах хх, (6) 
1=0 х; 


`Ця умова гарантуе ммум одне повне коливання зш(ах) на [а, 6]. 
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де 
К. хх <х,, 
ь УТ; + Цх- х, яет(АУ, ), Зах, 
Г <) = р х<х<х 1, 
1 > 
йе вх хм} 
Ум- м о 
И ЛЬ АУ; Е АУ; 
реек о: @ 


М- 


[№ 


Для похибки обчислення /(6) на клас! С г. м мае мсце ощнка знизу 
4Г 


г :(===) „г 
5Ш 5Ш х 
4[, 


оо ки 9 елы) [+26 6) 


[®/л] : 
2, дл 4 у ай «А! Е: УТ. 
@| л @+л @ 9 4 п 


«Погана функцщя» класу матиме вигляд: 


де Р(0) визначено в (5). 


Х + (хх; )ет($т @х,), хЕ Я 
55-м + Цх-хитим + 
+ (1+ зат (АЛ, (Л + Ех — хат (вт ох; Хер 


Ла + Их — хат (т ах; }, хе р, 


де хх, визначен! в (7). 


Тестування якост1 алгоритмив та ощнок 1х характеристик 
Для зручност! програмування формули (3) 1 (6) подамо у вигляд! 


М-1 
©) = Ул (оз р соз(ек;112 ) 
1=0 


М-—1 


В (®)= -(® соз(@ху ) + [9 - х; мет (м. соз(вх в 


1= 


|=) 


_ + [яет(АУ, {: —х; ое дров, } 
ГявтАУ;)[. [= 
+ т (ох, | эт (ох } — Ум соз(ах м )) 

@ 
ВДПОвВДНО. 

Тестування алгоритмав проводилося для сталих функщй, лиййних, кусково-лй- 
них, розглядалася «погана» задача класу Ст,м . Розглянемо деяк! з них. 

При тестуванн! обчислювалися точне значення /(@) та його наближення К(@), 
К («). Обчислювалися оцнки похибок (С, 62) 1 ИС; 02) для кламв Ср 1 Су ИДПОНДНО. 
Обчислювалися повн} похибки Е(©)= [(@)- Е(®) та Е! (@)= (®)- В! (®). 
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Приклад 1 
0, х<0,х>Ь 
2х, 0<х< 0,2, 
1 (х)=: -х+0,6, 0,2<х<0,4, хе 
0,2, 0,4<х<0,6, 
0,5х—0,1, 0,6<х<1 


Таблиця 1 — Результати теств для функцй Л (х) 


о М (со) ИСьо) Е (с) ИС мо) ЕК) 
2| 8001 800 | 2,26Е-04 | 3,90Е - 04 | 2,41Е- 04| 2,82Е — 04 | 113Е — 06 
2| 8001 1600 | 2,26Е-04 | 1,98Е- 04 | 1,17Е- 04| 1,44Е — 04 | 7,64Е — 07 
2| 8001 4000 | 2,26Е-04 | 7,98Е — 05 | 4,53Е - 05| 5,79Е — 05 | 5,63Е — 07 
2| 8001 8000 | 2,26Е-04 | 3,99Е — 05 | 2,22Е-05| 2,25Е — 05| 6,80Е — 04 
2| 8001 240001 2,26Е-04 | 1,33Е- 05 | 7,05Е-06| 3,38Е — 06 | 6,74Е — 04 
2 | 8001 80000 2,26Е-04 | 3,99Е - 06 | 1,80Е— 06| 2,90Е — 06 | 4 4ЛЕ- 07 
2 11000 1000 | -2,27Е- 04| 3,13Е- 04 | 1,30Е— 04| 2,27Е — 04 | 1,16Е — 04 
2 11000 2000 | -2,27Е- 04| 1,59Е-04 | 7,41Е-—05| 1,15Е — 04| 5ЛЕ- 07 
2 11000 5000 | -2,27Е- 04| 6,38Е — 05 | 311Е-05| 4,63Е — 05| 6,74Е — 07 
2 | 1000 10000] -2,27Е- 04| 3,19Е— 05 | 1,55Е —05| 2,32Е — 05 | 7,28Е — 07 
2 |1000 30000] -2,27Е 04 | 1,06Е-— 05 | 4,70Е —06| 7,7ЗЕ — 06 | 7,65Е — 07 
2 | 1000 100000 —2,27Е 04 | 3,19Е-06 | 8,68Е — 07| 2,32Е — 06 | 7,77Е — 07 
Приклад 2 
2 (х) =х зта(х) ‚ ХЕ [0,2 | 
Таблиця 2 — Результати тестування для функцй > (х) 
Е |® № (<) ИС) Е(с) ИС: м6) |Е\(@) 
2,5 1800 1800 |473Е- 04 7,40Е — 03 |4,/73Е 04 |7,40Е— 03 |2,54Е — 03 
2,5 800 11600 |4,/73Е — 04 199Е — 02 '4,/73Е-04 |1,03Е 02 |4,80Е — 04 
2,5 1800 14000 |4,73ЗЕ — 04 ЬПЕ- 02 14,/73Е- 04 |7,30Е— 03 |4,67Е — 04 
2,5 |800 18000 |4,73Е — 04 9,38Е —03 |4,73Е - 04 |7,25Е— 03 |4,69Е — 04 
2,5 800 124000 |4,73Е — 04 8,05Е — 03 14,73Е 04 |7,34Е— 03 14,72Е — 04 
2,5 |800 180000 |4,73Е — 04 7,60Е —03 |4,73Е 04 |7,38Е — 03 |4,73Е — 04 
2,5 1100011000 |-4,72Е 04 |743Е-— 03 14,72Е 04 |743ЗЕ— 03 \1,18Е - 03 
2,5 10002000 |-4,72Е- 04 |1,/74Е- 02 |4,72Е— 04 9/76Е— 03 14,68Е — 04 
2,5 10005000 |-4,72Е 04 |1,04Е- 02 |4,72Е— 04 |7,35Е— 03 14,76Е — 04 
2,5 100010000 |-4,72Е — 04 |9,01Е- 03 14,72Е — 04 7,31Е- 03 14,75Е - 04 
2,5 100030000 | -4,72Е 04 |7,95Е-— 03 |4,72Е 04 |7,38Е— 03 14,7ЗЕ -— 04 
2,5 11000100000 —4,72Е —04 |7,59Е — 03 |4,72Е 04 7АЛЕ- 03 14,72Е — 04 
Приклад 3 


Погана функц класу Сьм визначаеться значеннями {х;}, {1}, [. Як значення {7} 
взьмемо значення функцй /(х)= хзй(х). Результати тестйв подано в табл. 3. 
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Таблиця 3 — Результати тестування для функцй /(х)= хзш(х) 


8С 


Е |® № (<) ИС) Е) ИС, ме) Е! (©) 

2,5 |800 1800 4,73Е— 04  |7,40Е-— 03 |4,73Е — 04 |7,40Е — 03 |2,54Е — 03 
2,5 |800 |1600 [4,/73Е— 04  |1,99Е- 02 |4,73Е —04 |1,03Е — 02 |4,80Е — 04 
2,5 |800 |4000  |4,/73Е-—04  |1,ПЕ-02 |4,73Е — 04 |7,30Е — 03 |4,67Е — 04 
2,5 |800 |8000 [4,/73Е—04  |9,38Е— 03 |4,73Е —04 |7,25Е — 03 |4,69Е — 04 
2,5 |800 |24000 |4,/73Е—04  |8,05Е-— 03 |4,73Е —04 |7,34Е — 03 |4,72Е — 04 
2,5 |800 |80000 |4,/73Е—04  |7,60Е— 03 |4,73Е —04 |7,38Е — 03 |4,73Е — 04 
2,5 1100011000 |-4,72Е-04 |7,43Е— 03 |4,72Е —04 |7,43Е — 03 |1,18Е — 03 
2,5 1100012000 |-4,72Е-—04 |1,74Е— 02 |4,72Е — 04 |9,76Е — 03 |4,68Е — 04 
2,5 1100015000 |-4,72Е-—04 |1,04Е— 02 |4,72Е — 04 |7,35Е — 03 |4,76Е — 04 
2,5 11000110000 |-4,72Е-—04 |9,01Е-— 03 |4,72Е —04 |7,31Е — 03 |4,75Е — 04 
2,5 11000130000 |-4,72Е-—04 |7,95Е — 03 |4,72Е — 04 |7,38Е — 03 |4,73Е — 04 
2,5 |1000 1100000 |-4,72Е —04 |7,59Е — 03 |4,72Е — 04 |7,41Е — 03 |4,72Е — 04 


Вйдпов!дно до технологИ тестування [2] був проведений анал!з тестових таблиць. 
Результати тестування шдтверджують, що занурення задач! в Слыш вузький клас дае 
можливтсть б1лыш точно обчислити й розв’язок. Була шдтверджена яксть квадратурних 
формул та алгоритмав, що 1х реалзують, а також теоретичних ощнок похибки методу. 

За результатами тестувань можна зробити висновок, що точесть на практичних 
задачах вдрзняеться в середньому на порядок вд точности! на «поганйй» задач! класу. 

При збльшенн! числа вузлв вибтрки точнисть спочатку зростае, але з досягненням 
деякого № починае погршуватися за рахунок накопичення обчислювально! похибки. 
Несшвпадання апостерюрно! ощнки повно! похибки з ощнками (4), (8) пояснюеться вра- 
хуванням в апостерюрних ощнках похибки заокруглення. 
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ЕГ. Супроткин 

Тестирование качества вычисления оценки преобразования Фурье финитных функций 

В статье проведено тестирование качества алгоритмов вычисления преобразования Фурье на классах 
функций, соответствующих условию Липшица, и соответствующих им интерполяционных классах и 
продемонстрировано улучшение качества вычислений путем погружения задачи в более узкий класс. Были 
получены апостериорные оценки точности для рассмотренных тестовых задач. Подтверждено путем 
тестирования качество квадратурных формул и априорных оценок их точности. В статье приведены 
фрагменты таблиц с результатами тестов. 


Юотап биртойап 

Еошчег Тгап$Югт СасШаНоп Опайбу Тезбпо ое Ешйе Еипсвоп$ 

Тре {езипе оЁЕоипег гапзКопл сасшайоп а2огил$ даа ш Фе с1аззез сотезроп4ие ве ГлрзВи $ соп@ оп апа 
ш фе сотезропате пиегро!айоп с1аз3ез 15 ргезеще4; пиргоуетепе оЁ сасиайбоп диаШу Бу шса4те а ргоМет 
ш®ю а пато\уег с1а5$ 15 Четопзгае4. А розепоп езвтаноп оЁассигасу Юг ехатте рго Мега уаз оМашеа. Те 
ЧпаШу оР диайгаеаге Копли]а$ ап а рпоп езйтаез оЁ Фе ассигасу \еге сопйптед у {езйпе. ТВеге аге збоут 
Пастлет о ез соташие Фе {е5Ипе геи. 


Стаття надййила до редакций 21.07.2008. 
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